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We experimentally investigate the first-order correlation function of a trapped Fermi gas in the 
two-dimensional BEC-BCS crossover. We observe a transition to a low-temperature superfluid phase 
with algebraically decaying correlations. We show that the spatial coherence of the entire trapped 
system can be characterized by a single temperature-dependent exponent. We find the exponent at 
the transition to be constant over a wide range of interaction strengths across the crossover. This 
suggests that the phase transitions in both the bosonic regime and the strongly interacting crossover 
regime are of Berezinskii-Kosterlitz-Thouless type and lie within the same universality class. On the 
bosonic side of the crossover, our data are well-described by the quantum Monte Carlo calculations 
for a Bose gas. In contrast, in the strongly interacting regime, we observe a superfiuid phase which 
is significantly influenced by the fermionic nature of the constituent particles. 


Long-range coherence is the hallmark of superfluidity 
and Bose-Einstein condensation [1, 2]. The character of 
spatial coherence in a system and the properties of the 
corresponding phase transitions are fundamentally influ¬ 
enced by dimensionality. The two-dimensional case is 
particularly intriguing as for a homogeneous system, true 
long-range order cannot persist at any finite temperature 
due to the dominant role of phase fluctuations with large 
wavelengths [3-5]. Although this prevents Bose-Einstein 
condensation in 2D, a transition to a superfluid phase 
with quasi-long-range order can still occur, as pointed out 
by Berezinskii, Kosterlitz, and Thouless (BKT) [6-8]. A 
key prediction of this theory is the scale-invariant behav¬ 
ior of the first-order correlation function gi{r), which, in 
the low-temperature phase, decays algebraically accord¬ 
ing to gi{r) oc r~^ for large separations r. Importantly, 
the BKT theory for homogeneous systems predicts a uni¬ 
versal value of 77c = 1/4 at the critical temperature, ac¬ 
companied by a universal jump of the superfluid density 
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Several key signatures of BKT physics have been 
experimentally observed in a variety of systems such 
as exciton-polariton condensates [10], layered magnets 
[11, 12], liquid "‘He Aims [13], and trapped Bose gases 
[14-20]. Particularly in the context of superfluidity, the 
universal jump in the superfluid density was measured 
in thin Aims of liquid ^He [13]. More recently, in the 
pioneering interference experiment with a weakly inter¬ 
acting Bose gas [14], the emergence of quasi-long-range 
order and the proliferation of vortices were shown. 

There are still important aspects of superfluidity in 
two-dimensional systems that remain to be understood, 
which we aim to elucidate in this work with ultracold 
atoms. One question is whether the BKT phenomenology 
can also be extended to systems with nonuniform den¬ 


sity. Indeed, if the microscopic symmetries are the same, 
the general physical picture involving phase fluctuations 
should be valid also for inhomogeneous systems. How¬ 
ever, it is not known if algebraic order persists at all in 
the presence of inhomogeneity and particularly, whether 
the correlations in the whole system can still be charac¬ 
terized by a single exponent. Another fundamental issue 
that arises in the study of superfluidity is the pairing of 
fermions. While fermionic superfluidity has been exten¬ 
sively investigated in 3D systems [21-23], there are open 
experimental questions in the 2D context. In particular, 
what is the long-range behavior of spatial coherence of a 
2D fermionic superfluid, and can it also be described in 
the BKT framework like its bosonic counterpart? 

In this work, we probe the first-order correlation func¬ 
tion gi (r) of a trapped Fermi gas in the two-dimensional 
BEC-BGS crossover regime [24, 25]. The correlation 
function is determined from a measurement of the in situ 
momentum distribution of the gas. We demonstrate that 
even in this inhomogeneous system, algebraic order per¬ 
sists in gi(r) below a critical temperature. Furthermore, 
a quantitative analysis of the scaling exponents across 
the crossover reveals the validity of the BKT theory also 
in the fermionic regime. 

Our measurements are performed with a gas of 10® ®Li 
atoms confined in a highly anisotropic potential. The ax¬ 
ial and radial trapping frequencies are Wz « 27r x 5.5 kHz 
and Wj. « 27r X 18 Hz, leading to an aspect ratio of approx¬ 
imately 300:1. Our experimental system and methodol¬ 
ogy have been described in detail in Ref. [24] . We perform 
in situ imaging of the sample as a function of temperature 
and interaction strength. From the central density, we 
define the Fermi momentum fcp and Fermi temperature 
Tp, which constitute the relevant scales in the system. 
As shown in Ref. [24], for our experimental parameters. 
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all the relevant energy scales are smaller than the ax¬ 
ial confinement energy huiz- Hence the system is in the 
quasi-2D regime. 

We tune the interparticle interactions using a Fesh- 
bach resonance located at 832 G. Using the 3D scattering 
length asu [26] , the axial oscillator length £z [27] , and the 
Fermi momentum, we construct the effective 2D scatter¬ 
ing length a 2 D and crossover parameter ln(/cFa 2 D) [25]. 
For ln(fcFa 2 D) ^ — 1 and ln(A:Fa 2 D) 3> 1 we are in the 
bosonic and fermionic limit of the crossover, respectively. 

In addition to the measurements, we perform path- 
integral quantum Monte Carlo (QMC) computations of 
a Bose gas [28, 29] in a highly anisotropic 3D trap with 
parameters similar to those employed in the experiment. 
In the simulations, the bosons interact via the molecu¬ 
lar scattering length Omoi = 0.6 osd [30]. The relevant 
parameters that describe the system in terms of point¬ 
like bosons are the effective bosonic coupling strength 
g = y/^amo\/£z and the condensation temperature of an 
ideal 2D Bose gas « 140 nK, where 

TV is the number of particles. We use these bosonic 
parameters to compare our measurements to QMC at 
the lowest magnetic field values, where we have g = 
0.6,1.07,2.76,7.75 [31]. From the QMC computations, 
we obtain the local density prohle and the one-body den¬ 
sity matrix /9i(x,x') = (x)^(x')) for different inter¬ 

action strengths and temperatures, where (^(x) is the 
bosonic field operator. 

The global off-diagonal correlations in the system are 
encoded in the momentum distribution of particles. To 
reliably measure the in-plane momentum distribution 
h(k) of our sample, we employ the matterwave focusing 
technique described in Refs. [16, 32, 33], where the gas ex¬ 
pands freely in the axial direction while being focused by 
a harmonic potential in the radial plane. After expansion 
for a quarter of the period of the focusing potential, the 
initial momentum distribution is mapped to the spatial 
density profile, which we then image. We combine this 
focusing method with a rapid magnetic field ramp into 
the weakly interacting regime. This rapid ramp tech¬ 
nique - along with the fast axial expansion due to the 
large anisotropy of the trap - ensures that inter-particle 
collisions during the focusing do not cause significant dis¬ 
tortions to the measured momentum distribution. From 
h(k), we extract the absolute temperature T by means 
of a Boltzmann fit to the high-k thermal region [34]. 

To quantitatively investigate the spatial coherence in 
our system, we determine the first-order correlation func¬ 
tion gi (r) by means of a 2D Fourier transform of the mea¬ 
sured h(k). It is related to the one-body density matrix 
pi(x,x') by means of 



= J d^Rpi{R-r/2,R + r/2). (1) 



FIG. 1. First-order correlation function gi{r) for different 
temperatures at ln(fcFa 2 D) — —0.5 (upper left panel) and 
ln(fepa 2 D) — 0.5 (lower left panel). The temperature scale 
used here is t = T/T^^q. At high temperatures correla¬ 
tions decay exponentially as expected for a gas in the normal 
phase. At low temperatures, we observe algebraic correla¬ 
tions {gi(r) oc with a temperature-dependent scaling 

exponent g{T). This qualitative change of behavior is clearly 
visible in the for both exponential and algebraic fits (right 
panel), where a small value signals a good fit. In particular, 
this allows for an accurate determination of the transition 
temperature Tc (vertical dashed lines) [31]. 


A derivation of these relations is given in the supplemen¬ 
tal material [31]. The function (j'i(r) is a trap-averaged 
function, which captures the off-diagonal correlations of 
all particles in the system. Similarly, one can also define 
the central correlation function Gi(r, 0) = ((/)'!'(r)(^(0)), 
measured in the interference experiments [14, 35], which 
characterizes the correlations only in the central region of 
the trap, where the density is approximately uniform. In 
general, the two functions do not contain the same infor¬ 
mation and are only equivalent in a translation invariant 
system [31]. Note that due to the radial symmetry of the 
trapping and focusing potentials, the correlations only 
depend on distance and therefore it suffices to consider 
the azimuthally averaged function gi(r). 

Fig. 1 shows the experimentally determined gi{r) 
for different temperatures in the strongly interacting 
crossover regime. The correlation functions are normal¬ 
ized such that 5i(0) = 1. As expected, at high tempera¬ 
tures gi{r) decays exponentially with correlation lengths 
on the order of the thermal wavelength (At ~ 1.5/rm). 
As we lower the temperature, we eventually observe the 
onset of coherence over an extended spatial range that 
corresponds to several radial oscillator lengths £r, with 
£r ~ 6.8/rm. This shows that phase fluctuations in 
the system are non-local and span regions of the sam¬ 
ple where the density is not uniform. As pointed out in 
Refs. [36, 37], such extended spatial coherence in an in¬ 
teracting system is a sufficient condition for superfluidity 
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FIG. 2. Power-law scaling exponents across the two-dimensional BEC-BCS crossover. The temperature-dependent scaling 
exponent ri{T) in (a) the bosonic limit and (b) the crossover regime is shown. The relevant temperature scales in these cases 
are given by Tbec Tf, respectively. The crossover parameter ln(fcFa 2 D) is mildly temperature dependent. For reference we 
display the value at the critical temperature. For g = 0.60 (ln(fcFa 2 D) — —7.3) we show the prediction from QMC calculations 
for a Bose gas (filled red triangles) and an estimate of the effect of the hnite imaging resolution present in the measured data 
(open red triangles) [31]. We hnd an exponent which increases with temperature in agreement with BKT-theory. The power- 
law decay eventually ceases at Tc where a maximal exponent r/c is reached, (c) The value of r/c is approximately constant for 
all ln(feFa 2 D) where we have previously observed condensation of pairs [24]. This strongly suggests that the associated phase 
transition are within one universality class. 


in two-dimensional systems. 

As the temperature is lowered below a critical value, 
we find that the correlation function in an intermediate 
range 3At < r < 20At is well-described by a power- 
law decay, whereas exponential behavior is clearly disfa¬ 
vored. We quantify this by extracting the for both fit 
functions at different temperatures and observe a clear 
transition from exponential to algebraic decay (see Fig. 1 
b). This qualitative change in gi{r) provides an alterna¬ 
tive way to determine the phase transition temperature 
Tc from the kink in x^(T') [31]. We find that the cor¬ 
responding Tc obtained in this manner agrees with the 
temperature associated with the onset of pair condensa¬ 
tion that was measured in our previous work [24]. 

The power-law decay of gi{r) means that the spatial 
coherence of the entire sample is characterized by a sin¬ 
gle exponent 77 . Fig. 2 shows the experimentally deter¬ 
mined 77 for all the interaction strengths accessed in this 
work. We find g(T) to increase with temperature until it 
reaches a maximal value at T^, indicating a slower fall- 
off of correlations at lower temperatures. Although such 
temperature-dependence is qualitatively consistent with 
BKT theory, we observe the values of the exponents to 
be in the range 0.6 — 1.4 for the temperatures accessed 
in the measurement, which is substantially above the ex¬ 
pectation of ?7 < 0.25 for the homogeneous setup. 

To confirm the large scaling exponents in the trapped 
system, we compute the one-body density matrix on the 
bosonic side using the QMC technique described above. 
This allows to determine both the trap-averaged correla¬ 
tion function gi (r) as well as the central correlation func¬ 
tion Gi(r, 0). The trap-averaged gi{r) shows the same 
behavior as in the experimental case, i.e. a transition 


from exponential to algebraic decay at low temperatures. 
The corresponding QMC transition temperatures also 
agree with the measured values of Tc for g = 0.60,1.07, 
and 2.76. Furthermore, the maximal scaling exponent 
at Tc extracted from the QMC-( 7 i(r) for g = 0.6 is ap¬ 
proximately 1.35, which is close to the experimentally 
determined 7?(Tc) — 1-4. The central correlation function 
Gi(r, 0) shows a transition to algebraic order as well ~ 
with the same Tc as in the experiment - but with a max¬ 
imal exponent of approximately 0.25, as expected for a 
homogeneous system. This finding is also in agreement 
with the measurement of Gi (r, 0 ) in the interference ex¬ 
periments [14], and is explained by the nearly uniform 
density in the center of the trap. 

Fig. 2a shows the comparison between the experimen¬ 
tal and QMC values of g(T) for g = 0.60 (ln(fcFa 2 D) — 
—7.3). Although both show similar dependence on tem¬ 
perature, we find a considerable quantitative deviation 
between them. As discussed in the supplemental mate¬ 
rial [31], this discrepancy can mostly be attributed to the 
effect of the finite imaging resolution in the measurement 
of 7T,(k), which leads to an apparent broadening at low 
momenta and thus overestimates the value of 77 . We show 
an estimate of this temperature-dependent effect on the 
exponents (open red triangles) in Fig. 2a. There may be 
other effects in the experiment that contribute addition¬ 
ally to the deviation, such as higher order corrections to 
the determination of g from the fermionic scattering pa¬ 
rameters, and density-dependent inelastic loss processes. 

The experimental and simulated data raise the ques¬ 
tion why correlations in the trapped system decay with 
a larger scaling exponent than in the homogeneous case. 
To elucidate the role of inhomogeneity, we consider the 
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FIG. 3. Peak phase space density T>o = uqXt obtained 
from the central density no. The panels (a) and (b) show 
experimental and simulated data (for bosons) for the cou¬ 
pling strengths g = 2.76 and g = 7.75, respectively. The 
vertical dashed lines indicate the corresponding critical tem¬ 
peratures obtained from the measured onset of algebraic or¬ 
der. (a) We find excellent agreement between experiment 
and QMC for g = 2.76, providing evidence that we realize a 
strongly interacting 2D Bose gas. We verify the applicability 
of T>c = ln(380/g) [38] at this interaction strength (horizontal 
dashed line), (b) For the stronger coupling g = 7.75, however, 
we find the bosonic simulations to deviate from the measured 
results, indicating fermionic superfluidity. 

bosonic field operator given by ^{v) ~ \/p(j) exp(ii^(r)). 
In this representation, it is clear that one contribution 
to the decay of gi{r) in Eq. (1) comes from the spatial 
variation of the superfluid density p(r). Using a lo¬ 
cal density approximation and assuming the superfluid 
density to have a Thomas-Fermi profile, we estimate a 
contribution of approximately 0.3 — 0.4 [31] to the effec¬ 
tive exponent. Still, this fails to explain the large expo¬ 
nents observed in the experiment and the QMC simula¬ 
tions close to Tc. This suggests that the increase in the 
effective exponents is predominantly due to phase fluctu¬ 
ations in the inhomogeneous system, whose spectrum is 
modified by the discrete level-structure of the harmonic 
trapping potential and the Thomas-Fermi profile of the 
superfluid. This inference is further supported by calcu¬ 
lations of phase fluctuations in a trapped 2D Bose gas 
at low temperatures [39], which indicate a trap-induced 
increase of the effective exponent by up to a factor of 
three. 

Our measurements of gi{r) and ri{T) across the two- 
dimensional BEC-BCS crossover provide a unique op¬ 
portunity to study BKT physics even in the fermionic 
regime. Fig. 2 displays the measurement of the scaling 
exponent across the crossover. Remarkably, we find that 
- despite varying the scattering length by several or¬ 


ders of magnitude - the maximal scaling exponent rjc 
at the transition shows no dependence on the interac¬ 
tion strength (see Fig. 2c). We note that the actual 
value of 77c — 1.4 might depend on parameters specific 
to the experiment, such as particle number and trap¬ 
ping frequencies. Nevertheless, the fact that rjc remains 
constant across the BEC-BCS crossover unambiguously 
shows that the long-range properties at the transition 
are independent of inter-particle interactions. This is ev¬ 
idence that all the observed transitions for different in¬ 
teraction strengths lie in the same universality class. In 
particular, it shows that, even as we cross over to the 
fermionic side (ln(fcFa 2 D) > 0), the observed transitions 
are of BKT-type. 

We now turn to a quantitative investigation of local 
properties of the system. This allows to benchmark our 
measurements with (i) the QMC results for point-like 
bosons in the same quasi-2D trapping potential as real¬ 
ized in the experiment and (ii) QMC calculations of the 
homogenous 2D Bose gas [38, 40]. For this we investigate 
the phase space density (PSD) 

V = nXl. ( 2 ) 

Herein, n is the 2D density of atoms in a single hyperfine 
state and X^ = 2^11^/MkBT is the thermal wavelength 
of bosons with M being twice the fermion mass. Note 
that n coincides with the density of dimers in the bosonic 
limit. 

We first consider coupling strengths g = 0.60,1.07, 
and 2.76 on the bosonic side of the crossover. Fig. 3a 
shows the comparison between the experimentally mea¬ 
sured and QMC-computed values of the PSD in the trap- 
center for g = 2.76. We find excellent agreement be¬ 
tween the two data sets. In particular, at Tc, the central 
PSD for all three g are found to agree very well with 
Vc = ln(380/5) derived for a homogeneous 2D Bose gas 
with weak interactions (horizontal dashed line) [38, 40]. 
This shows that the onset of algebraic correlations in the 
trapped system coincides with the local PSD in the center 
of the trap crossing the critical value of the homogeneous 
system [28]. 

As we further increase ln(fcFa 2 D), the effective boson 
coupling strength g becomes very large. For g = 7.75 
(ln(fcFa 2 D) — 0.5), we find substantial deviations between 
the experimental and QMC data for the PSD at low tem¬ 
peratures (see Fig. 3b). Moreover, our QMC calculations 
show that the associated 2D Bose gas is in its normal 
phase for all temperatures accessed in the experiment. 
In contrast, the measurements show a clear superfluid 
phase transition at this interaction strength, as shown 
in Fig. 1 (lower panel). This provides evidence for the 
crossover to a superfluid phase whose properties are not 
captured by a description that assumes point-like dimers. 

Both experimental and simulated data in the bosonic 
limit are obtained in a highly anisotropic 3D trapping 
potential. Still, local observables such as the central 
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PSD and the central correlation function Gi (r, 0) agree 
excellently in their critical properties with the theory 
of a homogenous 2D Bose gas and the corresponding 
BKT phenomena. In the case of global correlations, we 
showed that the inhomogeneity leads to significant devi¬ 
ations from the homogeneous case, most importantly an 
increase in the exponent of the power-law decay. How¬ 
ever, the general features in the off-diagonal correlations 
- such as the temperature-dependence of ri{T) and the 
independence of rj^ from ln(fcpa 2 D) “ suggest that the 
long-range physics are still captured by the ideas under¬ 
lying BKT-theory for the two-dimensional XY model. 

In conclusion, we investigated the nature of the phase 
transition of a trapped 2D ultracold Fermi gas. We mea¬ 
sured for the first time the first-order correlation function 
of the entire system and extracted its long-range behav¬ 
ior. We showed that it is consistent with a description 
by a single power-law exponent for large distances. The 
transition temperature for onset of algebraic order coin¬ 
cides with the one obtained from the onset of pair con¬ 
densation in [24]. By comparing the experimental data 
to QMC calculations on the bosonic side, we found the 
system to realize a strongly interacting 2D Bose gas. The 
measured phase space densities and correlations on the 
fermionic side, instead, are not captured by a descrip¬ 
tion in terms of point-like bosons, which indicates the 
crossover to a fermionic superfluid. 

Our measurements show that the spatial coherence 
even in trapped systems can be characterized by a 
single scaling exponent. However, understanding the 
underlying mechanism remains a challenge for future 
explorations, and may lead to a deeper understanding of 
phase transitions in inhomogeneous systems. 
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SUPPLEMENTAL MATERIAL extracting critical temperature 

AND CRITICAL SCALING EXPONENT 


EXTRACTING THE FIRST-ORDER 
CORRELATION FUNCTION 


From the matterwave focussing technique we obtain 
the momentum distribution 

^(k) = (4«k) (3) 

of the trapped gas. Herein, dl is the creation operator 
for a particle with wave vector k. Note that h(k) involves 
an average over all particles in the trap. 

We now show that the Fourier transform of this func¬ 
tion coincides with gi(r) in Eq. (1). For this we recall 
that Ok is obtained from ^(x) by means of 

(4) 

Inserting this definition into (1) with the substitution 
s = R-|-r/2we arrive at 


ffiU = 


J (fs pi(s,r-f s) 

J d^s -f s)) 

d^fc (fk' 


d^s 


(27r)2 (27r)2 

(fk d^/c' 




(27r)2 (27r)2 


(27r)^(5^^^(k — k')e‘*' '^(akOk' 




(5) 


The expression in the last line is the Fourier transform 
of h(k). 

In a translation invariant situation, the one-body den¬ 
sity matrix can be written as pi(x, x') = /(x — x'), 
with some function /. In this case we have gi(r) = 
f d^s pi(0,r) oc Gi(0, r). The difference between the 
trap-averaged correlation function, 5 i(r), and the central 
correlation function, Gi(0,r), then consists of an overall 
(volume) factor, which vanishes in the normalization pro¬ 
cedure. 

Momentum resolution and coherence length: 

The matterwave focusing lens used for the measure¬ 
ment of ?T.(k) has a magnification factor Mwiensj where 
‘k'lens = 27r X 10 Hz is the trap frequency of the focus¬ 
ing potential and M is the molecular mass. From this, 
the effective momentum space resolution can be obtained 
according to Ak = McuiensAx/h ~ 0.035^m“^, with 
Ax ^ 5 /im being the spatial imaging resolution. This 
means that the largest coherence length that is accessi¬ 
ble in gi(r) is approximately L = 271 /Ak ^ 105 fj-m. 


The qualitative change of the decay of correlations 
when lowering the temperature allows for a determina¬ 
tion of the critical temperature for each value of a 2 D- 
For this purpose we fit both an exponential {gi{r) = 
06 “'"/'’) and algebraic {gi{r) = ar~^) model function to 
the intermediate length scales of gi (r) and extract the as¬ 
sociated x^-value. A smaller value of corresponds to a 
better fit. In Fig. S 1 (upper panel) we display X^(T) for 
812G {g = 7.75). A sharp transition in the behavior of 
correlations is visible at a certain temperature, which we 
associate with the critical temperature Tc- To determine 
Tc we piecewise linearly interpolate X^{T) according to 


Xaig(U = ci0(Tii) -T) + C 2 {T - Tii))0(r - 

( 6 ) 

xLp(U = C3(ri2) - T)0(Ti2) -T) + c^e{T - Ti^)) (7) 

for the algebraic (alg) and exponential (exp) fits, respec¬ 
tively. Here 9{x) is the Heaviside step function. We 
generically find ci and C 4 to be small, which justifies the 
choice of a power law ht at low temperatures, and an 
exponential fit at large temperatures. Furthermore, the 
temperatures and coincide within a few percent. 



Fig. S 1. Extracting critical temperature Tc and critical ex¬ 
ponent gc at 812 G. The upper panel shows the y^-values for 
exponential (red) and power law (blue) fits of gi{r), respec¬ 
tively. Lower values of indicate a better description of the 
data. This allows to determine the critical temperature as the 
temperature of onset of algebraic decay in gi{r). The lower 
panel shows the scaling exponent g{T) below Tc. We deter¬ 
mine the critical exponent (red circle) by extrapolating g{T) 
using a polynomial fit. 

















B [Gauss] 

g 

\n{kFa2D)Ta 

(stat.)(sys.) 

^c/^bec (stat.)(sys.) 
(algebraic decay) 

Tc/Tf (stat.)(sys.) 
(algebraic decay) 

Tc/Tf (stat.)(sys.) 
(pair condensation [24]) 

Vc (sys.) 

692 

0.59 

-7.30 (4) (1^) 

0.695 (25) (Igf) 

0.091(21) {t\l) 

0.089 (15) (lit) 

1.35 (12) 

732 

1.07 

-3.42 (2) (1^) 

0.590 (24) (l)f ) 

0.094(19) (1)^) 

0.100 (22) (lit) 

1.36 (18) 

782 

2.76 

-0.59 (1) (1"^) 

0.483 (20) (Ifg) 

0.114 (19) (t?^) 

0.129 (35) (HI) 

1.48 (25) 

812 

7.75 

0.57 (1) (1?) 

0.485 (22) (IPs) 

0.142 (24) (t?|) 

0.146 (25) (li) 

1.42 (12) 

832 

- 

1-23 (1) (1^) 

- 

0.157 (30) {tlD 

0.167 (39) (1«) 

1.33 (18) 

852 

- 

1.72 (1) (1") 

- 

0.166 (40) (tfa) 

0.167 (27) (til) 

1.50 (15) 


TABLE 1. We show the data for the measured critical temperatures and critical exponents for all interaction strengths considered 
in this work. For each magnetic field value we display the corresponding bosonic coupling strength g and the crossover parameter 
ln(fepa 2 D) at the critical temperature. The Fermi momentum is defined as fep = \/inn, where n is the central density of atoms 
in a single hyperfine state. The critical temperature obtained from the onset of algebraic decay is shown with respect to Tbec 
and Tp. For comparison, we also list the transition temperatures for the onset of pair condensation determined in [24]. The 
scaling exponent 77 c = v{Tc) is obtained from the extrapolation described in Sec. II. 


We set Tc = (T^'^ +ri^^)/2 to obtain the transition tem¬ 
perature. 

Using the critical temperature Tc found in this man¬ 
ner, we extract the scaling exponent rjc at the transition. 
For this purpose we extrapolate the experimental data 
points for r](T) by means of a quadratic polynomial fit, 
see Fig. S 1 (lower panel) and extract the value of r](Tc). 
We list the measured critical temperatures and scaling 
exponents in Table I. 

Errors: The method described above to extract the 
critical exponents contains some uncertainties. The sta¬ 
tistical errors on the measured exponents are quite small 
and the error on rjc as shown in Fig. 2.c is mainly due 
to the uncertainty in the extrapolation of rj{T). The 
error bars are obtained according to 5rjc = \r](Tc^'^) — 

SYSTEMATIC EFFECTS 

The errors on experimental quantities shown in the 
main text are statistical uncertainties of our measure¬ 
ments. The systematic effects in our measurements have 
been discussed in detail in [24] . In this work, we addition¬ 
ally introduce the temperature scale T/Tg^Q and extract 
the scaling exponent r]{T) from gi{r). These quantities 
are systematically affected by the following factors: 

Uncertainty in particle number 

We determine the number of atoms in the cloud from 
in-situ absorption images. N is affected by the inten¬ 
sity of the imaging beam, magnification of the imaging 
system and the small population of atoms in the non¬ 
central pancakes of the trapping potential. The ideal gas 
condensation temperature depends on the atom number 


according to = ^6N , and hence it is affected 

by the uncertainty in JV. In addition to experimental un¬ 
certainties, the measurements also contain atom number 
fluctuations of about 10-15%. 

Imaging effects 

The correlation function gi(r) is determined by means 
of a 2D Fourier transform of the in-plane pair momen¬ 
tum distribution h(k). The momentum distribution is 
obtained using a matterwave focusing technique which 
consists of a ballistic expansion of the gas in a harmonic 
potential for a quarter of the trap period (t/4 = 25 ms) 
and subsequently imaging the planar density distribu¬ 
tion. 



Fig. S 2. Effect of finite imaging resolution on the scaling 
exponent g for g = 0.60. The imaging simulations are per¬ 
formed assuming the 31 (r) and the corresponding exponents 
g obtained from QMC computations (filled triangles). The 
open triangles show the exponent extracted after the simu¬ 
lated imaging. We find a substantial temperature-dependent 
deviation that is qualitatively consistent with the experimen¬ 
tally measured exponents (filled circles). 
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The onset of algebraic decay at large distances in gi{r) 
corresponds to a peak at low momenta in the measured 
h(k). Naturally, a peakier h(k) leads to a broader gi{r) 
and hence a smaller scaling exponent. This method of 
determining gi{r) and rj{T) from the momentum distri¬ 
bution is limited by two main factors: 

a. Vertical expansion during TOF: As shown in 
[33], the gas expands rapidly in the vertical direction 
upon release from the trapping potential. After 25 ms, 
the vertical extent of the cloud is approximately 500 /im. 
This can lead to some parts of the sample exceeding 
the depth of focus of the imaging system, which in turn 
causes some distortion in the measured density distribu¬ 
tion n(r,T/4). In general, it leads to a broadening of 
the fi(k) which has the effect of increasing the measured 
scaling exponent 77 . 

b. Finite imaging resolution: As we lower the 
temperature, the momentum distribution becomes nar¬ 
rower. The measured h(k) is the convolution of the ac¬ 
tual momentum distribution with the hnite resolution of 
our imaging setup. This convolution leads to a broaden¬ 
ing of the momentum distribution for small k and hence 
a steeper decay of gi{r) for large r. Moreover, this broad¬ 
ening effect is larger for distributions that are closer in 
width to the resolution limit. Intuitively, this means that 
the distortion caused by the imaging resolution is en¬ 
hanced at lower temperatures. For the extracted scaling 
exponents, this leads to a temperature-dependent devi¬ 
ation between the actual value and the measured value. 
The scaling exponent is always overestimated due to the 
finite imaging resolution. 

To estimate the contributions of the vertical extent of 
the sample and the finite imaging resolution on the ex¬ 
tracted scaling exponents, we perform a simulation of the 
imaging setup which consists of two lenses (/ = 80 mm) 
in a 4f-conhguration and a CCD camera. 

From the QMC-computed gi{r) at each temperature, 
we determine h(k) by means of an inverse Fourier trans¬ 
form. This corresponds to a 2D spatial distribution 
fi(r) of the cloud with jrj = Ii|k|/Mwiens) where wiens 
is the trap frequency of the focusing potential. To ac¬ 
count for the vertical extent of the cloud after time-of- 
flight, we construct a 3D density distribution according 
to n{x, y, z) = L~^ ■ h(x, y) for —L. 2,12 < z < Tz/ 2 , where 
Lz « 500 ym is the vertical size of the cloud. While this 
is not truly reflective of the actual distribution in the ex¬ 
periment, it is sufficient to capture the essential effect of 
the vertical size of the gas. 

We consider the imaging effect of a thin section of the 
cloud defocused by a distance 2 ; from the focal plane at 
z = 0. Using the paraxial wave equations, we compute 
the propagation of the imaging beam, approximated as a 
plane wave, through this section and through the lenses. 
The resulting intensity distribution contains the effect of 
the defocus as well as the finite resolution of the lenses. 
We perform this computation for all —Lz/2 < z < Lz/2 


and integrate the resulting intensity distributions in the 

L /2 

4f plane according to hf{x,y) = ^J^j^I{x,y,z)dz. This 
yields the imaged column density n'{x,y) and hence the 
imaged momentum distribution h'(k). From this we 
obtain the imaged g'i{r) and the corresponding scaling 
exponent Fig. S2 shows the comparison of ex¬ 

ponents extracted from the experimental measurement, 
QMC computations and the imaging simulations, for typ¬ 
ical simulation parameters. 

From the imaging simulations, it is clear that the finite 
imaging resolution causes a signihcant overestimation of 
the scaling exponent. We find that the main contribu¬ 
tion to the deviation is in fact from the finite resolution in 
the radial plane and the effect of the vertical extension of 
the cloud is mild. This explains the strong temperature- 
dependence of the discrepancy between measured and 
QMC exponents as shown in Fig. S 2. 


LOCAL DENSITY APPROXIMATION 

We estimate the influence of the Thomas-Fermi (TF) 
profile of the superfluid density on the correlations within 
a local density approximation. The result of the analy¬ 
sis can be summarized in the following statements: (1) 
The TF profile results in an additive contribution ?7 tf 
to the scaling exponent, which is approximately 0.3 at 
low temperatures; (2) The temperature-dependence of 
ryxF is mild and from the QMC data for the density pro¬ 
files we estimate the value of ?7 tf at the transition to be 
approximately 0.4. We conclude that the large scaling 
exponents observed in the experiment, in particular the 
value rjc — 1-4 at the transition, are mostly due to phase 
fluctuations in the inhomogeneous sample. 

We assume that, within a phase-amplitude represen¬ 
tation, the bosonic held in the superhuid phase can be 
written as 

<^(r) = \4^expii^(r), (8) 

where 4> is an operator, but p{v) is a function. We ap¬ 
proximate p{v) to be given by a TF prohle according to 
p(v) ~ Po(l ~ — ''’^/.Rtf); where 6{x) is the 

Heaviside step function, and Atf = {2ti^gno)^^'^/Mujr is 
the radius of the superhuid core. The particular shape 
of p(r), however, is not essential for the following conclu¬ 
sion. We then hnd 

5 i(r) = J p{r)p{s -H (9) 

for the trap-averaged correlation function. 

Approximating the phase huctuations to be translation 
invariant we write ~ |r — with 

a temperature-dependent exponent ?7phase(F), which is 
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assumed to be constant throughout the superfluid region. 
We arrive at 

( 10 ) 

itTF 

with the function / being shown in Fig. S 3a. 

The spatial decay of the function / results in an 
additive contribution Srjes = —dln//dln(r/i?TF) — 
1.97(r/i?TF)^ to the measured scaling exponent. In our 
experiment, we have i?TF — 100 /rm-^/g and typical fit¬ 
ting ranges are r < rgt — 4£r = 30 /rm. Defining the TF 
contribution as ?7tf = Rtf) we have 

V{T) = rjTFiT) + VphnseiT) ( 11 ) 

for the total scaling exponent extracted from the data. 
The temperature dependence of the TF contribution re¬ 
sults from the temperature dependence of the central 
density no(T). We use the latter from the QMC data for 
the density profiles to compute ??tf(T") = ^tf (To) j 
where Tq is a reference temperature. We choose the lat¬ 
ter to be small. The estimated TF contribution 77tf(T) 
for g = 0.60 is shown in Fig. S 3b. 



Fig. S 3. Estimated effect of the Thomas-Fermi (TF) profile 
onto the decay of correlations in a local density approxima¬ 
tion. Panel (a) shows the function f{r/ Rtf) from Eq. (10) 
which multiplies the algebraic decay due to phase fluctuations. 
Our fitting range, highlighted by the blue shaded region, is 
given by r / Rtf 0.4, where the function is rather fiat. The 
fall-off of /(f/Ttf) leads to an additive contribution to 
the extracted scaling exponent. The latter is shown in panel 
(b). At low temperatures, where phase fluctuations are small, 
we have g ~ t;tf — 0.3. Even at large temperatures, we only 
have tjtf — 0.4. This effect cannot fully explain the large ex¬ 
ponents found from the QMC and experimental data. Hence 
the associated decay of correlations must be mostly due to 
phase fluctuations. 











